Macroscopic quantum fluctuations in noise-sustained optical patterns 
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We investigate quantum effects in pattern-formation for a degenerate optical parametric oscillator 
with walk-off. This device has a convective regime in which macroscopic patterns are both initiated 
and sustained by quantum noise. Familiar methods based on linearization about a pseudo-classical 
field fail in this regime and new approaches are required. We employ a method in which the pump 
field is treated as a c-number variable but is driven by the c-number representation of the quantum 
sub-harmonic signal field. This allows us to include the effects of the fluctuations in the signal on 
the pump, which in turn act back on the signal. We find that the non-classical effects, in the form 
of squeezing, survive just above the threshold of the convective regime. Further above threshold the 
macroscopic quantum noise suppresses these effects. 

PACS number(s): 42.50.Lc, 42.50.Dv, 42.50.Ct, 42.65.Sf 



I. INTRODUCTION 

Nonlinear optics has provided an ideal testing ground 
for ideas in both nonlinear dynamics and quantum op- 
tics. It provides fast nonlinearities and a degree of con- 
trol that allow fundamental dynamical systems to be re- 
alized and nonlinear phenomena like pattern formation 
to be demonstrated [jj. It also provides systems with 
very low levels of noise so that fluctuations can be lim- 
ited by quantum effects. The combination of these fea- 
tures has led to the study of quantum phenomena in op- 
tical pattern formation and of noisy precursors of 
the patterns, which have been termed quantum images 
||-|3| . The accurate modeling of such quantum nonlinear 
systems presents a significant challenge. Pattern forma- 
tion and dynamics is usually associated with excitation 
of a large number of transverse modes and a fully quan- 
tum description of each of these is required in order to 
properly treat the quantum fluctuations. The Heisen- 
berg picture produces a hierarchy of coupled nonlinear 
operator equations that usually defies analysis. The pre- 
ferred method to date has been to linearize the quantum 
fluctuations about a classical field amplitude that usually 
takes a constant value below threshold HQ , but may be 
associated with a stable pattern above threshold Q . 

A more difficult situation arises if the system displays 
macroscopic features driven by noise. In such cases we 
cannot expect linearization of the quantum fluctuations 
to give reliable results and a new approach is needed. A 
simple device demonstrating macroscopic, noise-driven 
patterns is the degenerate optical parametric oscillator 
in the presence of walk-off. The semi-classical analysis 
of this device reveals a region of convective instability, 
above the threshold for oscillation, in which noise sus- 
tained structures are seen in the transverse field distri- 
bution The aim of this paper is to develop a suit- 



able approximation scheme with which to model quan- 
tum effects in the parametric oscillator in this regime of 
operation. 

The convectively unstable regime is characterized by 
an amplification and flow of fluctuations [fll[ . In systems 
in which the spatial reflection invariance is broken by the 
presence of a group velocity term, local perturbations of 
the steady state can be advected more rapidly than their 
growth rate. If the system is deterministic then at any 
fixed point any initial localized perturbation decays and 
the system approaches the undisturbed steady state. In 
this case macroscopic patterns can arise and be observed 
only if noise is continuously applied, the structure now 
being regenerated at any time, hence the name noise- 
sustained patterns. These structures are the result of 
noise self-organization, with magnification factors of sev- 
eral order of magnitudes. They are thus interesting can- 
didates for the study of quantum correlations in spatially 
structured systems. 

Any system with an advection (or drift or walk-off) 
term that is also not translationally invariant will, in 
general, be convectively unstable when operating suffi- 
ciently close to and above the onset of the instability of 
the steady state. Hence, this type of instability has been 
predicted in a number of optical systems including Kerr 
media with a tilted pump [p|Jr2]| and optical parametric 
oscillators (OPO) with walk-off 

Modeling quantum effects in the regime of convective 
instability for a nonlinear optical device presents a double 
challenge. First, the system has a broad spectrum both 
in frequency (at a fixed point) and in wave- vectors (far 
field at a fixed time), thus it cannot be studied within 
a few-mode approximation. Second, we should be able 
to follow the evolution of the fluctuations from the mi- 
croscopic level through the amplification into the macro- 
scopic pattern. In order to do this we introduce, in Sect. 
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IV, a suitable non-linear approximation with which to 
treat the convective regime of a degenerate optical para- 
metric oscillator. In order to fix the terms of reference for 
this approximation we begin, in Sect. |n|, with a review of 
the semi-classical features of the device and its convective 
instability. This is followed, in Sect, [Iljj, by a quantum 
description of the device. Once we have introduced our 
method we discuss quantum features of the device in its 
various regimes of operation (Sect. |vj), paying particular 
attention to the demanding convective regime (Sect. VI). 



II. SEMI-CLASSICAL DESCRIPTION AND 
CONVECTIVE REGIME 

We consider a Degenerate Optical Parametric Oscilla- 
tor (DOPO), a device consisting of a cavity filled with 
a nonlinear medium, which converts a pump at fre- 
quency 2lo into a sub- harmonic signal at frequency uj. 
The possibility of phase matching the down-conversion 
process depends on the birefringence of the crystal that 
provides a difference of refractive index for differently po- 
larized fields. We can exploit this difference in order to 
avoid the effects of dispersion by selecting the same in- 
dex of refraction for the pump and signal n 2LJ = n u . In 
this paper we consider type-I phase matching for which 
ordinary polarized pump photons are down-converted to 
produce pairs of extraordinary polarized photons that are 
degenerate both in frequency and in polarization. 

In anisotropic media rays do not necessarily travel in 
a direction perpendicular to their wavefronts pE[ . As a 
consequence the extraordinary-polarized signal generated 
in our DOPO will walk off, that is it will propagate in the 
transverse direction relative to the ordinarily-polarized 
pump. This transverse walk-off effect is described in the 
dynamical equations by a term which accounts for a ve- 
locity relative to the frame of reference fixed by a pump 
of finite transverse width. 

The quantum effects we wish to study are associated 
with the convective regime and it is important to define 
this carefully. The different regimes of operation of a 
DOPO can be understood within a semiclassical theory 
and this section provides a brief (semi-classical) analysis 
of the convective and other regimes. A more complete 
discussion can be found in p3| . The intra-cavity field 
is described by two slowly varying complex field ampli- 
tudes Ao(x,t) and A\(x, t) for the pump and signal re- 
spectively. These depend on the transverse spatial coor- 
dinates x — (x, y) and the time t. Within the paraxial 
approximation (for propagation in the z direction), the 
mean field limit and for single longitudinal mode opera- 
tion the dynamical equations become 
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d t A {x,t) 
d t Ax(x,t) 



- 7o [l + iA - ia W 2 ]A (x, t) - 9 -A\{x, t) 

+E (x) + e £ (f, t), 
-71 [1 +£Ai -iaiV 2 -vd y ]Ai(x, t) 



(1) 



Here £j (i = 0, 1) are additive Gaussian white sources of 
noise, with non-vanishing correlations of the form: 



< t)£f (f, >= dijdix - x' )6(t - t'). 



(3) 



The level of noise introduced is fixed by the parameters 
eo and t\. Our fully quantum analysis will produce equa- 
tions of similar form in which these parameters are fixed. 
Eq is the amplitude of the driving field which we take 
to be real. The remaining parameters in these equations 
are the cavity decay rates 7,, the cavity detunings Aj, the 
diffraction aj, the walk-off v and the nonlinear coefficient 
g. It is convenient to introduce scaled variables 



t' = -ft 



(4) 



E' 



7 3/2 fl D/4 



where we have restricted the cavity decay rates and 
diffraction coefficients such that 7 = 70 = 71 and 
a = ao = ai/2. Our equations are valid either for one or 
two transverse spatial dimensions (D — 1,2). On omit- 
ting the primes, our amplitude Eqs. (Q) and (^J) become 

d t A (x, t) = -[l + iA - iV 2 }A a {x , t) - -A\{x, t) 

+E (x) + e $o(x,t), (5) 
d t Ai (x, t) = -[l + iA 1 - 2iV 2 - vd y ]A 1 (x, t) 

+A {x,t)A* 1 {x,t) + e 1 ^ 1 (x,t). (6) 

For a uniform driving field Eo, the equations (|^) and (||) 
admit the homogeneous stationary solution 



An 



E () 



1 + iAo 



A\ = 



(7) 



The threshold for parametric oscillation can be deter- 
mined by a linear stability analysis of this solution. The 
linearized equations for signal and pump fluctuations 
5Ai(x,t) — Ai(x.t) — Af (i = 0,1) are decoupled, and 
the fluctuations of the pump are always damped. For the 
signal, we consider perturbations of the form e?k-£+\{k)t 
and find the dispersion relation 



A±(fc) = -l + ivk y ± V^ 2 - (Ai +2|fc| 2 ) 2 , (8) 



where we have introduced a scaled pump 

Eo 



F = 



a 2 , 



(9) 



We find that there is an instability at F = 1. For 
F < 1, Re(X) < and the solution (0) is absolutely 
stable. For F > 1, there is a positive growth rate of 
fluctuations (i?e(A+) > 0) which takes a maximum value 



+gA Q (x 1 t)A\{x, t) + ei&(2, t). 



(2) for \k c \ = yj— Ai/2 if the signal detuning is negative 
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(Ai < 0), and for k = if Ai > 0. In this paper we 
are interested in the case of pattern formation and we 
restrict our analysis to the case Ai < 0. The instability 
at F = 1 when v — is a Turing instability in which a 
stationary pattern appears |L6| . If v ^ then the eigen- 
value becomes complex and we find a Hopf bifurcation in 
which a traveling pattern emerges Jl3| . 

The direction of instability is determined by the eigen- 
functions V±(k, —k) of the linear problem dtV±(k, —k) = 
X±(k)V±(k, —k). Solving this gives 



V±(k,-k) = e i * ± 5A 1 (k)±5A$(-k) 

iAi + 2z|fc| 2 T v / l^l 2 -(Ai + 2|fc| 2 ) 2 



(10) 



e i4> ± (fe) _ 



The solution V-\-{k, — k) gives the direction of amplifica- 
tion of fluctuations, while fluctuations are damped for 
V—(k,—k). In particular, for the critical wave-vector 
\k c \ and for a real pump Eq and Ao = 0, we obtain 
V±(k c ,-k c ) = 5Ai{k c ) ± 5A\{-k c ). Therefore, in this 
case, the difference of real parts and the sum of imagi- 
nary parts of field in k c and — k c will show damped fluc- 
tuations at threshold. We also note that the instability 
direction is independent of the walk-off term. 
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FIG. 1. Stability diagram as a function of the signal de- 
tuning Ai. The different lines correspond to the threshold 
of absolute instability (F c ) for different values of the walk-off 
parameter: v — 0.2 (dotted line), v — 0.42 (continuous line), 
v = 0.6 (dashed line). When F < 1 the solution (Q) is abso- 
lutely stable, while for 1 < F < F c the solution is convectively 
unstable. 

Above the instability threshold (F=l) the steady state 
(Q) is convectively unstable: any perturbation grows 
while travelling in the direction fixed by the walk-off term 
and eventually leaves the system |]l3| . In this regime a 
continuous perturbation such as a source of noise, gives 
rise to a noise-sustained pattern consisting in disordered 
travelling stripes in the signal. On increasing the pump 
a second threshold is reached at F = F c . Beyond this 
threshold the pattern is sustained by the nonlinear dy- 
namics, being also present in the absence of perturba- 
tions, once it is formed. The state (R) is absolutely un- 
stable in this regime |13fl. In Fig. [if we plot the result 



of the calculation of the absolute instability threshold as 
function of the signal detuning Ai, for different values of 
the walk-off parameter v. 

Walk-off has three main effects in this process of pat- 
tern formation jl2|,[l3]]. The first is the existence of the 
convective regime in which patterns are sustained by 
noise. Second is that it breaks the rotational symme- 
try, favoring the formation of stripes orthogonal to the 
walk-off direction and traveling in this direction. Thirdly 
the selected wave-vector, that is the most intense mode 
A:m of the pattern, depends on the walk-off parameter. 
An approximate expression for kM c & n be obtained in 
the context of front propagation into an unstable state 

There are two important characteristics of the noise 
sustained patterns that exist in the convective regime. 
The first is a broad spectrum, both in frequency and in 
wave- vectors [ [l9| . Second is the presence of macroscopic 
amplified signal fluctuations around the unstable ref- 
erence state (0). These characteristics imply that the 
convective regime cannot be studied within a few-mode 
approximation, because many modes contribute signifi- 
cantly to the spectral properties. The presence of macro- 
scopic fluctuations also invalidates approximations based 
on linearization schemes. These facts make a quantum 
formulation of the convective regime especially difficult. 
We face a situation in which nonlinearities determine the 
dynamics of fluctuations around the reference state, with 
fundamental quantum noise being amplified by several 
orders of magnitude to produce a macroscopic pattern in 
the signal. 



III. QUANTUM FORMULATION 

In the quantum formulation of the DOPO the intracav- 
ity pump and signal fields are given by operators A (x, t) 
and A\ (x, t) that satisfy standard equal-time commuta- 
tion relations [pfl 



Ai(&,t),At(x",t) 



6ijS(x 



(11) 



where the indices i,j stand for 0, 1. Following the tech- 
niques described in || , we can introduce a model Hamil- 
tonian for the device. This will include the effects of 
diffraction together with the driving by a real, classi- 
cal external field, nonlinear interaction between the fields 
and cavity damping. Our model, however, also requires 
that we take account of the effects of walk-off. The re- 
sulting Hamiltonian gives, on making the usual Markov 
approximation, the coupled Heisenberg equations: 



d t A (x,t) = - 7o [(l + iA ) - m V 2 ]i (x,i) 
-|i 2 (f,t)+ E (x) + F a 

dtA^x, t) = - 7l [(l + iAi) ~ «*iV 2 - vdy}!^, t) 



+gA {x, t)A\(x, t) + Fi 



(12) 

) 

(13) 
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Note that these are very similar in form to the semi- 
classical Eqs. (||) and (g). The Langevin operators Ft 
describe the quantum noise added as a consequence of 
the interaction with the bath of external modes. These 
have the non-vanishing second moments: 

(F l (x,t)F J \x',t')) = 2 ll S lJ S(x-x')5(t-t'). (14) 

A direct solution of these non-linear Langevin equa- 
tions of operators is impractical, requiring the solution of 
an infinite hierarchy of equations for the evolution of all 
the products of operators that are coupled by the dynam- 



ics. A standard alternative approach to this Heisenberg 
picture is to consider the evolution equation of the re- 
duced density operator p of the system in the Schrodinger 
picture and to use quasi-probability functionals. In this 
approach to the quantum dynamics of open problems, the 
intracavity dynamics is described by a master equation 



(15) 



where H is the Hamiltonian 



H 



.3=0,1 



k-it(f)(Aj - a j V 2 )A ] {x) + ijxvAli^dyAxix) + iE (x)(Al(x) - A (x)) 



i y -{A {x)A\ (x)-h.c). 



The Liouvillian A accounts for dissipation through the 
partially reflecting mirrors of the cavity and is given by 



d 2 x^j 



{[A0)r P A\W] + [A0)p,A]{x)]). 



The master equation (|l5| ) can be mapped onto an equa- 
tion of motion for one of a number of quasi-probability 
distributions in the phase-space of the system [^o] p2[ . 
These distributions are functionals of the c-number fields 
ai(x) associated with the operators Ai(x). This evolu- 
tion equation is obtained by substituting products of field 
operators and the density operator, depending on the or- 
dering, by suitable operations on the distribution func- 
tionals i|il< 

The evolution equations obtained in this way for the 
distributions are functional partial differential equations. 
These are not in general of the Fokker-Planck type and 
do not lead to well-behaved stochastic representations in 
terms of Langevin equations driven by Gaussian white 
noise. In particular the Hamiltonian term describing the 
X 2 interaction gives a contribution: 



[A (x)A\ (x)~h.c.,p}< 



1 



SafScfQ 



Sa 



2aoa*- h c.c. 

oai 



where the parameter s depends on the ordering. This 
term does not fulfill the requirements that guarantee a 
positive definite solution for W s : in the Wigner repre- 
sentation (s = 0) we find third order derivatives, while it 
is known [p3| that positiveness requires a Fokker-Planck 
form of the master equation (only first and second order 
derivatives) or to include derivatives to all orders. For 
the P (s = 1) and Q (s — —1) representation third order 
derivatives disappear, but the diffusion matrix is not pos- 
itive definite so that positive solutions are again not guar- 
anteed, although the Q retains positivity through having 



a minimum allowed width [ pl| . Generally these prob- 
lems have been avoided by using linearization schemes 
p2j. Such linearization approximations, however, are 
valid only for small damped fluctuations. They cannot 
be used in a convective regime as the reference state is 
unstable and the fluctuations, far from being small, are 
amplified. The alternative of the P positive represen- 
tation |24|] is not suitable for the same reason and the 
unstable reference state results in diverging trajectories. 

These problems of the convective regime can be illus- 
trated by comparison with the situation of a DOPO be- 
low the threshold of signal generation. In this case the 
stable solution is a homogeneous pump with an ampli- 
tude that depends on the coherent driving field. The 
signal field is zero on average, but its fluctuations show a 
level of self-organization that increases near the thresh- 
old. This is the regime of quantum images noisy 
precursors generated by quantum noise. These images 
reflect the presence of eigenmodes of the linearized equa- 
tions, whose eigenvalues are such that their negative real 
part approaches zero at threshold. The fluctuations of 
these eigenmodes are the least damped ones and domi- 
nate the dynamics of the signal. The important point is 
that the intensity of such quantum images of the signal 
is of the order of quantum noise, while the pump has a 
macroscopic mean value. It is then possible to neglect the 
fluctuations in the pump, approximating it by a classical 
coherent field p^ j. In this approximation the Hamilto- 
nian is a quadratic function of the quantum operators. 
The consequence is that a well defined Fokker-Planck 
equation for the Wigner distribution is obtained. Such 
Fokker-Planck equations can be represented in terms 
of stochastic Langevin equations for the c-number field 
cti(x) |6|. The same type of approximation, linearizing 
around a pattern solution || , is generally possible in the 
absolutely unstable regime above threshold. A common 
feature of these two regimes (absolutely stable and unsta- 
ble) is that quantum noise does not change drastically the 
solution with respect to the stable deterministic solution. 
This means that in the stochastic representation, fluctu- 
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ations only induce the trajectory to visit a small region in 
phase space in the neighborhood of the deterministic so- 
lution. In the convective regime the classical determinis- 
tic solution is unstable and macroscopically different from 
the stochastic solution. In this regime quantum noise in 
the DOPO is amplified, destroying the zero-valued ho- 
mogeneous deterministic solution for the down-converted 
field and driving the system into noise sustained states 
having a macroscopic number of photons. 

IV. TIME DEPENDENT PARAMETRIC 
APPROXIMATION 

In this section we we propose an approximate descrip- 
tion of the quantum dynamics of the DOPO in the con- 
vective regime, based on the main physical features of 



this regime. Our aim is to be able to treat the macro- 
scopic quantum fluctuations associated with the signal 
field in the convective regime. 

In the convective regime there are large signal fluctu- 
ations around the unstable solution A\ = 0. The cou- 
pling of signal and pump gives the nonlinear saturation 
for these amplified fluctuations. On the other hand, the 
pump field is always macroscopic and stable, with small 
damped fluctuations. This suggests the approximation of 
neglecting quantum noise in the pump and approximat- 
ing it by a classical field Aq (x, t) . In this way we obtain 
a Hamiltonian that is quadratic in the operators describ- 
ing the quantum dynamics of the signal field. For such 
quadratic Hamiltonians, the Wigner quasi-probability 
functional of the complex function ot\{x, t) obeys the 
following Fokker-Planck equation in which the classical 
pump field Ao(x,t) appears parametrically: 



S \ S 2 

— 71 [(1 -HAi) -miV 2 - vd y ]ai(x,t) + gA (x,t)al(x,t) + c.c.J + 7i 5ai5a * W(ai;A )). 



dW(a i; Ao) 
dt 

The associated Langevin equation that represents the 
stochastic dynamics of the signal field ot\ [x, t) is 

d t ax(x,t) = -7i[(l + iAx) - iaiV 2 - vd v ]ai(x, t) + 

gAo(x, t)a\(x,t) + ^/ji£i(x,t), (16) 

where £i (x, t) is a complex Gaussian white noise (see 
Eq.(|])). This noise term accurately represents the ef- 
fects of vacuum fluctuations associated with cavity losses 
on the signal field. We note that treating the pump field 
classically in this way is a natural extension of the para- 
metric approximation to three-mode interactions, which 
treats a strong mode classically and has been widely used 
in quantum optics for many years Q . 

It is important to note that Aq(x, t) cannot be replaced 
by an expectation value of (Aq) as would be possible in 
the regime of absolute stability (quantum images). Such 
an ansatz decorrelates the pump modes from the sub- 
harmonic ones and eliminates the saturation effect of the 
pump. In fact, with such an ansatz Eq. ( |l6|) becomes 
linear, giving a Gaussian probability distribution for the 
signal modes. This distribution would always be cen- 
tered on zero, but with statistical moments that diverge 
above threshold because the signal modes are undamped 
in the convective regime. Therefore, the stochastic dif- 
ferential equation must be solved self-consistently with 
an equation defining the dynamics of the classical field 
Aq. The equation we propose for Ao is suggested by the 
Heisenberg equation (|lj) , with Aq replaced by a classical 
field Aq. We first neglect the noise source in (JIJ) since 
quantum fluctuations entering in the cavity are unim- 
portant, as compared with the macroscopic fluctuations 
of the signal term A\. Secondly we replace the oper- 
ator A 2 by the c-number function a\ associated with 
our stochastic representation of the signal. This replace- 
ment is independent of operator ordering and hence will 



be the same should be choose to use a different quasi- 
probability. This procedure gives a partial differential 
equation for the "classical" pump field driven by the c- 
number representation of the quantum signal field: 

d t A (x, t) = -7o [(1 + iA ) - ia W 2 }A (x, t) 

-|a 2 (f,t) + E (x) (17) 

A justification for this equation is that its mean value co- 
incides with the expectation value for Aq obtained from 
the operator equation (p"2|). This procedure is reminis- 
cent of the time-dependent refinement of the parametric 
approximation described in some detail by Kumar and 
Mehta |2(| . This approach allows for the quantum evolu- 
tion of the weak fields to feed back and affect the classical 
strong field. In the approach of Kumar and Mehta, this 
feedback is via quantum expectation values of operators 
for the weak fields. Here, however, we are required to take 
explicit account of the noisy properties of the quantum 
sub-harmonic field. We do this by using the c-number 
representation of the quantum field, associated with our 
stochastic simulation of it, as a term in equation (|l7j). 

In summary, our time dependent parametric approxi- 
mation is defined by stochastic classical equations in the 
Wigner representation for the fields Aq and ot\, which, 
with the scaling (^), are: 

d t AQ{x,t) = -[{l + iA )-tV 2 }Ao(x,t) 

-±al(x,t) + E a (x) (18) 
d t ai{x,t) = - [(1 + iAi) - 2iV 2 - vd y ] ai(x,t) 

+A Q (x,t)a* 1 (x,t) + -^ 7I ^Z 1 (x,t). (19) 

Stochastic averages of the c-number variable a\(x,t) will 
provide symmetrically ordered averages of the quantum 
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fluctuations in the signal field as driven by the "classi- 
cal" pump field. The classical pump field is driven by the 
macroscopic quantum fluctuations in the signal as rep- 
resented by the c-number representation a%(x, t). This 
time dependent parametric approximation appears use- 
ful in situations in which there are large fluctuations of 
the signal that cannot be described by approximations 
based on linearization. 



V. STOCHASTIC TRAJECTORIES AND 
WIGNER DISTRIBUTION FUNCTION 

Numerical simulation of the stochastic trajectories as- 
sociated with the Langevin equations © and © gives 
a good intuitive understanding of the dynamical proper- 
ties of the regime below threshold, the convective regime 
and the absolutely unstable regime. In this section we 



present such numerical simulations working with a single 
transverse dimension (D — 1) ^7j. 

Fig. H is a space-time plot of the near-field for the 
signal in the below-threshold, convective and absolutely 
unstable regimes. Fig. ^ gives the far-fields associated 
with the same simulations. In the following we discuss 
the different properties of these trajectories and how they 
are reflected in the associated Wigner distribution. In 
particular, we consider the phase space dynamics of the 
most intense modes of the signal pattern. The Wigner 
probability distribution associated with these modes dis- 
plays distinctly non-Gaussian features in the convective 
regime. These are a result of the interplay of non-linear 
and walk-off effects. It is clear that they cannot be de- 
scribed within a linearization scheme that does not take 
account of this interplay between the signal and pump 
fields. 



X 






FIG. 2. Evolution of the near field of the real part of the signal Re(ai(x, t)) for: (a) F = 0.999, (6) F = 1.025, (c) F = 1.1, 
in 2500 time units. Parameters are: Aq = 0, Ai = — 0.25, « = 0.42, system size = 1.7678 • 512 pixels ~ 900 space units. Only 
the regions in which the signal is excited are shown, that is the region of the plateau of the supergaussian pump fegj. 



T I 



k 



I. 



(6) 



(o) 



FIG 3. Evolution of the far field \a x (k, t)\ for: (a) F = 0.999, (b) F = 1.025, (c) F = 1.1. Same parameters as in Fig. | 



A. Below threshold 

Below but close to threshold we find weakly damped 
fluctuations which are a precursor to the traveling pat- 
tern that appears at threshold. The fluctuations impose 
a degree of spatial self-organization in those regions in 
which the pump is sufficiently strong to bring the OPO 
close to threshold. In Fig. ||(a) we plot the real part 
of the stochastic variable ai(x,i) for a single trajectory. 
This is a realization of these fluctuations for a pump with 
a supergaussian profile |28| . Noisy patterns of this form 
have been predicted for the bclow-threshold OPO with- 
out walk-off and have been termed quantum images . 
Not too close to threshold, the damped fluctuations can 
be analyzed with linearization procedures Q in the limit 
of small fluctuations. Our nonlinear quantum equations 
enable us to study also the regime closer to threshold, 
where large critical fluctuations are expected to occur. 
Note, in particular, that the results in Fig. ||(a) were 
obtained for F = 0.999. 

The selection of a preferred wave-number in the 
stochastic pattern of Fig. Wa) becomes more evident 
in the far field shown in Fig. ||(a). It is clear that there 
are preferred values of the wave-number but that a broad 
distribution of weakly damped modes around these pre- 
ferred k c and —k c modes is apparent. 



space trajectory for the slowly varying amplitude of the 
dominating Fourier component, a'i(k c , t), is shown in Fig. 
|(a). 

The linear stability analysis of Sect. || also identified 
the direction of instability V+. In particular, in the case 
of a real pump, and for the critical mode k ci this direction 
is given by [ai(+fc c )+a 1 (— k c )\. As a consequence the su- 
perposition of modes [a.\{+k c ) + k c )\ can be decom- 
posed in two quadratures, one corresponding to the di- 
rection of instability that becomes undamped at thresh- 
old (Re[ai(+k c )+ai(— k c )]), and the orthogonal quadra- 
ture (Im[cti(+k c ) + ai(— k c )]) that remains damped. We 
observe that the superposition of slowly varying modes 
a[(±k c ,t) can be decomposed into damped and un- 
damped quadrature in the same way. In fact due to 
the symmetry u){k) = — u;(— k) we have V±(k, — k) — 
e i "W t [e i *±5A' 1 (lc)±5A'?(-k)], so that the relative phase 
e l *± between the slowly varying modes is the same as 
that in the equation (|i"o|). Hence, we can also identify 
the real and imaginary quadratures of the superposi- 
tion of modes [ai(+k c ) + a[(— k c )] as damped and un- 
damped at threshold. The corresponding time trajectory 
of this superposition of modes displays very clearly the 
expected reduction of fluctuations in the damped imagi- 
nary quadrature (see Fig. |](b)). 
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0.0 
Re[a\(kJ] 

FIG. 4. (a) trajectory of the slowly varying ampli- 
tude ai(+k c ) during 20.000 time units. (6) trajec- 
tory of [a' 1 (+k c ) + a'i( — k c )] during 20.000 time units. 
F = 0.999, other parameters are as in Fig. ^, except for 
dx = 51A C /512 ~ 1.7702, where A c = 2ir/k c , 512 is the num- 
ber of grid points. 

An interesting characterization of the stochastic dy- 
namics in the far field, Fig. ||(a), is obtained by looking 
at the time evolution of the stochastic amplitudes for the 
most intense modes a\(k c , t). We first recall that the lin- 
ear stability analysis of Sect. [n] identifies the existence 
of a non-vanishing frequency (w(k) = vk) at threshold 
caused by the walk-off. This implies that a traveling 
pattern will emerge above threshold and that the corre- 
sponding Fourier modes will oscillate at this frequency. 
We can remove this time-dependence by working in a 
frame rotating at this frequency. This corresponds to 
factoring out a time factor e 1 "*-^* to obtaining the slowly 
varying amplitude a[(k,t) — a\(k, tje - ^"'. A phase 




FIG. 5. Wigner distribution for the superposition of 
modes Qi(+fe c ) + k c ). Parameters of Fig.W. Total time 
2.000.000 units. Note the factor 10 difference in the scale of 
the two axes. 

From the stochastic trajectories that randomly visit 
the different points of phase space it is easy to construct 
a relative histogram giving a probability density in this 
phase space. This density is identified with the Wigner 
distribution. As with all Wigner functions, the marginal 
distributions, obtained for one field quadrature by inte- 
grating over the orthogonal quadrature, are true prob- 
ability distributions for the remaining quadrature. At 
a finite distance from threshold the Wigner distribution 
W{ot\{k)) for the field a\(k) obtained in this way has 
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a Gaussian shape consistent with a linearized analysis 
of fluctuations. Such a Gaussian Wigner distribution is 
a solution of the Fokker- Planck equation for the Wigner 
representation of linear signal fluctuations. If we consider 
the Wigner distribution for the superposition of modes 
discussed above W(ai(+k c ) + oti(—k c )), then we obtain a 
Gaussian centered on the origin but with a variance that 
depends on the orientation in phase space [29| . There 
is an axis with a reduced variance ('squeezed') and the 
orthogonal one with a larger variance ('anti-squeezed') 
(see Fig. ||). These features reflect the asymmetry or 
phase-sensitivity of the fluctuations already visualized in 
the stochastic trajectory. 

B. Convective regime 

Differences between the regime below threshold and 
the convective regime are clearly seen both in the near 
and far signal fields. We observe a macroscopic traveling 
pattern in the near field (Fig. ||(fo)). This is clearly asso- 
ciated with wave-numbers distributed around the value 
of the selected one (k M ) in the far field (Fig. The 
spectrum of excited wave-numbers is clearly narrower in 
the convective regime than below threshold. This is re- 
flected in the more regular pattern appearing in the near 
field. Our simulations display the typical features asso- 
ciated with the convective regime Jl3| : 

• The noise sustained pattern does not fill the whole 
region in which the pump has a value above thresh- 
old. This is because the pattern grows while trav- 
eling in the walk-off direction. Note that the space 
point at which the pattern reaches a macroscopic 
observable value changes randomly from time to 
time. This reflects the origin of the pattern in 
(quantum) noise. 

• The far field shows the predominance of different 
wave-numbers at different times resulting in a spa- 
tial spectrum that is broader than that found in 
the absence of walk-off or in the absolutely un- 
stable regime. There is competition between the 
modes within this broad spectrum and hence it 
is not possible to define, in this regime, a single 
wave-number Um corresponding to the most ex- 
cited modes. Modes with different wave-numbers 
compete to form the pattern, switching on and off 
as the pattern evolves. 

Phase space trajectories for this regime are shown in 
Fig. ^. We find that there are random changes in the 
phase and amplitude of the slowly varying signal a'^+kc) 
around a zero mean value (Fig. ||(a)). This is similar to 
the behavior depicted in Fig. [|(a) below threshold. The 
difference is that in the convective regime macroscopic 
intensities are reached, with the signal amplitude taking 
values comparable to those reached in the absolutely un- 
stable regime (compare scales of Figs. ||(a), ^|(a), and||). 



The continuous changes in intensity from zero to macro- 
scopic values originate in the fact that, in the convec- 
tive regime, a given mode is not constantly switched-on 
(see FigJ||(fe)). The pattern is sustained by noise and is 
subject to a continuous renovation: different stripe pat- 
terns (with different wave-numbers) grow, travel in the 
system starting from noise and die out. This has an im- 
portant consequence in the time scales of the far field 
dynamics: below threshold these scales are determined 
by noise, while in the convective regime they are deter- 
mined by the time needed for a perturbation to travel 
through the system. Another indication of the nonlin- 
ear dynamics of fluctuations that occur in this regime is 
that the quadrature displaying reduced fluctuations is no 
longer the one determined by the linear analysis. This is 
seen in Fig. ||(&) where the ellipse of fluctuations is tilted 
with respect to the corresponding one below threshold 
Fig. |(6). 




-50 50-100 100 

Re[a\(kJ] Re[a\(k c )+a\(-kJ] 



FIG. 6. (a) Trajectory of slowly varying ampli- 

tude of a'i(+k c ) during 100.000 time units. (6) Trajec- 
tory of [a' 1 (+k c ) + a[(-kc)]. Parameters F = 1.025, 
A = 0, Ai = -0.25, v = 0.42, dx ~ 1.7702. 




FIG. 7. W(ai(+fc) + ai(— k)), for an excited mode 
k = 1.04fc c , obtained from a trajectory during 10.000.000 time 
units. Other parameters as in Fig. pi 

The probability distributions obtained from the tra- 
jectories of Fig. ^ also reflect the nonlinear nature 
of the fluctuations in this regime. In Fig. |?] we 
show the W distribution for the superposition of modes 
[ai(+k) + ot\{— k)] for one of the most excited wave- 
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numbers, namely k — 1.04fc c . A most noticeable feature 
is the non-Gaussian shape of the distribution for large 
values of the amplitude in the direction of undamped 
fluctuations. The wings of the distribution originate in 
the macroscopic fluctuations of the mode under consider- 
ation when it switches-on. Its most probable value is still 
zero, reflecting the fact that most of the time the mode 
remains switched-off. We can view these nongaussian fea- 
tures in the wings of our Wigner functions as precursors 
of the pair of peaks appearing in the absolutely unsta- 
ble regime. These wings become more pronounced as we 
approach the absolutely unstable regime. 

Finally, we note that the modes that become excited 
and contribute to the dynamics seem to reach a com- 
mon maximum amplitude. This is probably fixed by the 
maximum value of the energy exchanged with the pump 
mode in the nonlinear interaction. This is shown in Fig.^ 
where the possible values of different modes are seen to 
be cut-off at essentially the same amplitude. The non- 
gaussian form of these distributions is also clear and this 
again demonstrates that we are dealing with nonlinear 
effects associated with the quantum fluctuations. 




-50 



50 



FIG. 8. Section of the Wigner distribution along the imag- 
inary axis W(Q, Im(ai(k))) for 3 excited modes: k c (dashed 
line), k' = 1.04fc c (dash dot line), k" = 1.06fc c (continuous 
line). Same parameters as Fig. 



tensity are maintained. Although there is essentially only 
the phase diffusion for &m, the critical mode, with wave- 
number k c , displays a second frequency superimposed on 
the phase diffusion process. 
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FIG. 9. Trajectories of (a) a'i(k c ) and (b) a^km) during 
100.000 time units, F = 1.05, other parameters as in Fig. ^| 
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FIG. 10. 



-200 200 

Re[a\(k u ) + a\(-k„)] 
Trajectory of [a'i(-l-fawr) + ai(— fawr)] over 



10.000.000 time units. Other parameters as in Fig|j. 



The phase space trajectory for the superposition of 
modes [a[(kM,t) + a^(— kM,t)] is shown in Fig. We 
observe that fluctuations are not uniformly distributed 
around a zero value as they were in the below threshold 
(Fig. [|) and convective (Fig. ^]) regimes. Instead, they 
describe a closed curve around the origin. The associ- 
ated W distributions display peaks at two values. These 
correspond to the two points of maximum curvature of 
the elliptical ring. 



C. Absolutely unstable regime 

In the absolutely unstable regime we observe from the 
near field plot, Fig. ||(c), that a macroscopic and stable 
traveling pattern fills the whole of the above threshold 
region. This behavior is reflected in the far field, Fig. 
|§|(c) , which shows a well-defined and fixed dominant wave 
number and a narrow spatial spectrum. We should note 
that the dominant wave-number /cm does not coincide 
with the most unstable wave- number at threshold (fc c ). 
This is a consequence of the interplay between nonlinear- 
ities and walk-off. Phase space trajectories for the am- 
plitudes of these two modes are shown in Fig. [)| Even 
after elimination of the rapid frequency there remains a 
phase diffusion process, but macroscopic values of the in- 
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FIG. 11. W(Re{ai(k)), 0), for positive values of Re[ai(k)], 
for (a): k = k c , (6): k = kM, obtained from a trajectory 
during 10.000.000 time units. Other parameters as in Fig. H. 
These figures are symmetric around 0. 
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The main characteristics of the trajectories in phase 
space are reflected in the associated Wigner distributions. 
For the less intense modes contributing to the dynam- 
ics we can approximate the associated Wigner function 
W{a\{k)) by a Gaussian, displaced from and orbiting 
about the origin in phase space. In Fig.|ll|(a) we show 
a cut along the real direction of the Wigner distribution 
for the critical mode (W(Re[ax(k c )], 0)). By contrast, 
the most intense mode (with wave- number &m) displays 
some interesting new features. Fig.|ll|(6) shows an asym- 
metry in the distribution of fluctuations around the mean 
amplitude in each of the peaks, with a sharp decay of the 
distribution at some maximum amplitude. These facts 
indicate the existence of nonlinear properties associated 
with the quantum fluctuations in the absolutely unstable 
regime. These features would necessarily be absent in 
any analysis based on a linearization about a determin- 
istic macroscopic state. 



VI. NON-CLASSICAL PROPERTIES IN THE 
CONVECTIVE REGIME. 



The convective regime is characterized by amplified 
fluctuations and macroscopic noisy patterns. It is inter- 
esting to ask, therefore, if any of the low-noise quantum 
features found below threshold can survive in this noisy 
environment. Quantum effects in the OPO have been 
observed as sub-shot noise fluctuations both in the field 
quadratures and intensity differences associated with the 
down-converted light j3(]] . Examples of the noisy features 
associated with the real part of the signal field in this 
regime are plotted in Fig. [l^ for three different values of 
the driving field, all within the convective regime. Note 
the different scales on the vertical axes in these figures. 
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FIG. 12. Snapshots of the real part signal Re(ot\(x)) 
for different pump values:(a) F = 1.001, (6) F = 1.01, (c) 
F = 1.025. Other parameters are: dx = 1.7678, 512 grid 
points, a = 1, A = 0, Ai = -0.25, v = 0.42. Note the 
different vertical scales in the figures. 

It is helpful, in looking for non-classical effects, to keep 
in mind the manner in which such effects appear below 
threshold. We will also restrict ourselves to the study 
of quantum correlations in the far field. Conditions for 
squeezing and associated non-classical effects are usually 
expressed in terms of normally ordered moments of op- 
erators (indicated by : :). These can be obtained from 
the symmetrically ordered moments (indicated by £()), 
that are associated with the Wigner function, by use of 
the commutation relations dill): 



(: A(k,t)A(k',t) :) = (S(A(k,t)A{k',t))) 

(: A(k,t)A\k',t) :) = (5(i(fc,<)i t (fc',i))) - ~5(k - k') 

The S function appearing in the second of these equations 
is a signature of the shot or vacuum noise. Our approxi- 
mation scheme is based in the Wigner representation and 
gives results for correlations of symmetrically ordered op- 
erators for the intracavity fields. In order to obtain re- 
sults for the corresponding normally ordered products 
and to test for the presence of non-classical effects, we 
need to establish a reference shot noise level. This level 
can be obtained for each quadrature correlation from the 
variance of the linear stochastic process associated with 
the empty cavity: 

d t s{x, t) = - [(l + iA!)-2zV 2 ] s{x,t) + -^ u ^ x {x,t). 

Here we have omitted the walk-off term as it does not 
affect the shot noise level. Squeezing in our simulations 
will be associated with a quadrature probability distribu- 
tion that is narrower than the Gaussian associated with 
this linear process. In general we can consider a different 
quadrature for each wave-number k. It is useful to define 
a pair of (superposition mode) quadratures for each k pa- 
rameterized by the angle 9. For the critical wave-number 
these take the form 



Ax{k c ,t) ± A x {-k c ,t) 



h.c 



(20) 
squeezed 



o 



300 



600 



900 



We expect, in general, that the most strongly 
quadrature should depend on the value of 8 [Q. 

We begin our investigation of the convective regime 
at a point that is just above threshold with F = 1.001 
(Fig. [l2](a)). Fluctuations associated with the pattern 
are in this case still relatively small and we find that the 
Wigner distribution has a Gaussian shape as shown in 
Fig.[l3|. We find that there is quadrature squeezing, with 
the squeezed quadrature X- (0) exhibiting the same level 
of squeezing as is found just below threshold. In partic- 
ular, for F — F t hr ± 0.001 we find that the intracavity 
field is squeezed by 50% below the shot noise limit for 
a flat pump and by 37% for a supergaussian pump [pl[. 



fO 



This indicates a smooth variation across threshold for the 
squeezed quadrature variance. For excited modes, other 
than the critical one, we also find squeezing below the 
shot noise level for the appropriate quadrature. 




-0.01 0.01 

^[a^kj-a^-kj] 

FIG. 13. W{Re{ai(k c ,t)-ai(-k c ,t)),0). Continuous line 
is obtained for F = 1.001, and dotted line for F = 1.01. The 
dashed line represents the distribution for the vacuum state, 
corresponding to the shot noise level. The distributions are 
relative to trajectories of 2.000.000 time units. Other param- 
eters as in Fig. 

Increasing the value of the pump, so as to move further 
into the convective regime, leads to a rapid increase in 
the magnitude of the signal field. Indeed, for (F = 1.01) 
we observe, in Fig. |l2](6) that the signal field has grown 
by two orders of magnitude. Fluctuations are still ex- 
tremely phase sensitive and, as depicted in Fig. ^(a) 
there is a strong reduction in the fluctuations for some 
quadratures. This reduction is insufficient, however, to 
reach below the shot noise level and there is no squeezing. 
In fact, we find residual fluctuations +27% above the shot 
noise level. This is comparable with the value associated 
with the coherent states. These enhanced fluctuations 
are associated with a much broader Wigner distribution 
as shown in Fig. [l^. It is remarkable, however, that this 
enhanced but still small level of fluctuation can coexist 
with the macroscopic fluctuations in orthogonal quadra- 
ture. If we move still further above threshold then we 
find, for F = 1.025 (Fig. |l2](c)), a variance which is 159 
times the shot noise level and both quadratures display 
fluctuations that are well above the level usually associ- 
ated with quantum effects. We note that for the parame- 
ter values used here, the threshold of absolute instability 
for an infinite system occurs at F ~ 1.035. 
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FIG. 14. Var(X_(0)) at F — 1.025 for the critical wave- 
number k — k c . The minimum occurs for 9 < 0. 



A further indication of the nonlinear nature of the fluc- 
tuations in the convective regime is given by the fact that 
the angle 9, for which there is the greatest reduction in 
the fluctuations, changes with the strength of the pump 
value. This has already been discussed in connection with 
Fig. |(a). In particular, for the critical wave-number, 
X_(9) shows strongest squeezing for 6 = in the lin- 
ear regime below threshold. In the convective regime, 
however, the greatest reduction in the quadrature fluc- 
tuations occurs for a value of < 0. This is shown in 
Fig. |lj in which we plot the variance V&r(X-(9)) in nor- 
mal ordering and normalized to the shot noise level, for 
F = 1.025. 

The OPO can also exhibit strong correlations be- 
tween the far field intensities associated with oppo- 
site wave-numbers. We have calculated the fluctu- 
ations in the intensity difference for opposite wave- 
numbers associated with the normally ordered moment 

(: (Ai(k)Ai(k) - Ai (-k)Ai(-kj) :)■ A negative va- 
lue for this quantity indicates a non-classical effect some- 
times referred to as twin beams or intensity-difference 
squeezing |32| . As in our discussion of quadrature squeez- 
ing, we find that this quantity is only negative very near 
to threshold (F = 1.001). Further into the convective 
regime we find that the macroscopic noise associated with 
the formation of a pattern increases the noise in the inten- 
sity difference. For F = 1.01 we find that the intensity- 
difference squeezing has been replaced by fluctuations in 
excess of the shot noise level. 

In summary, we have shown that quantum effects can 
survive above threshold in the convective regime but only 
very near to threshold. On increasing the pump and en- 
tering further into the convective regime, we find that 
non-linear effects associated with the fluctuations tend to 
distribute part of the macroscopic fluctuations into the 
observables that are squeezed nearer to threshold. This 
identifies walk-off as an effective mechanism of quantum 
decoherence in which the macroscopic nonlinear fluctua- 
tions present in the convective regime overwhelm quan- 
tum effects associated with noise reduction. 



VII. CONCLUSION 

We have introduced a suitable method to describe the 
quantum properties of macroscopic patterns sustained 
by quantum fluctuations in a degenerate optical para- 
metric oscillator with walk-off. These patterns appear in 
the convective regime and are characterized by a broad 
far field spectrum with continuous competition between 
several wave-numbers (thus, few mode approximations 
are not adequate) and by being the result of amplified 
quantum fluctuations around a unstable reference state. 
Traditional linearization techniques cannot be applied in 
these situations. Instead we use a time dependent para- 
metric approximation in which the pump field is treated 
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as a c-number variable but driven by the c-number repre- 
sentation of the quantum sub-harmonic signal field. The 
key point is that this includes the effects of the fluctua- 
tions in the signal on the pump which in turn act back 
on the signal. 

Using this method we have described the quantum fluc- 
tuations in type-I OPO with walk-off in three regimes: 
below the threshold of instability, in the convective un- 
stable regime and in the absolute unstable regime. 

Below threshold we find that the Wigner representa- 
tion has a Gaussian shape centered at the origin. This is 
the result previously found in a OPO without walk-off ||] 
from a linearized analysis. We also find that the walk-off 
does not destroy the existence of squeezing in suitable 
quadratures. 

In the convective regime the macroscopic character of 
the fluctuations is reflected in a extremely broad Wigner 
distribution where the probability is still centered at the 
origin but the nonlinear effects lead to the appearance of 
wings in the distribution which is no longer a Gaussian. 
These wings are in fact precursors of the pair of peaks 
appearing in the absolutely convective regime. We show 
that squeezing in the appropriate observables can be also 
obtained in this regime but only just above threshold. 
The walk-off and the nonlinearities act as quantum deco- 
herence mechanism, distributing part of the macroscopic 
fluctuations into the observables that were squeezed be- 
low threshold. Another nonlinear effect appears in the 
selection of the quadrature displaying reduced fluctua- 
tions, that is no longer the one determined linearly. 

In the absolutely unstable regime there are also clear 
indications of nonlinear properties associated with quan- 
tum fluctuations. The interplay between walk-off and 
non-linearity results in a complex dynamics in which the 
frequencies of far-field modes are not constant, giving a 
complicated variation of the phases. Also, the most in- 
tense mode is not the critical mode. We find that while 
the Wigner distribution for the less intense modes can be 
approximated by a Gaussian (displaced from the origin 
and orbiting about it) this is not the case for the most 
intense modes for which the distribution of fluctuations 
is asymmetric around the mean amplitude with a sharp 
decay at some maximum amplitude. 

Finally, our method can be used in other situations and 
systems in which there are large fluctuations of the signal 
that cannot be described by approximations based on lin- 
earization. This includes situations in which the critical 
fluctuations appear at threshold for pattern generation. 



ACKNOWLEDGMENTS 

This work was supported by the Quantum Structures 
Network of the EU TMR programme (Project No. ERB 
FMRX-CT96-0077). RZ, MSM and PC acknowledge fi- 
nancial support from MCyT (Spain) projects PB97-0141- 
C02-02 and BFM 2001-1108. SMB thanks the Royal So- 



ciety of Edinburgh and the Scottish Executive Education 
and Lifelong Learning Department for financial support. 



[1] S. Ducci, N. Trips, A. Maitre, and C. Fabre, Phys. Rev. 

A 64, 023803 (2001) and references therein. 
[2] L. A. Lugiato and F. Castelli, Phys. Rev. Lett. 68, 3284 

(1992). 

[3] L. A. Lugiato, A. Gatti, and H. Wiedemann, in Quantum 
Fluctuations, edited by S. Reynaud, E. Giacobino, and J. 
Zinn Justin (Elsevier-North-Holland, Amsterdam, 1997) . 

[4] A. Gatti and L. A. Lugiato, Phys. Rev. A 52, 1675 
(1995); L. A. Lugiato, A. Gatti, H. Ritsch, I. Marzoli, 
and G. L. Oppo, J. Mod. Opt. 44, 1899 (1997); I. Mar- 
zoli, A. Gatti, and L. A. Lugiato, Phys. Rev. Lett. 78, 
2092 (1997). 

[5] For a review, see for example (a)L. A. Lugiato, Chaos 
Solitons Fractals 4, 1251 (1994), and references quoted 
therein;(b) L. A. Lugiato, M. Brambilla, and A. Gatti, in 
Advances in Atomic Molecular and Optical Physics, 40, 
229, B. Bederson and H. Walther eds. (Academic Press, 
New York, 1998). 
[6] A. Gatti, H. Wiedemann, L. A. Lugiato, I. Marzoli, G. L. 

Oppo, and S. M. Barnett, Phys. Rev. A 56, 877 (1997). 
[7] J. Jeffers and G-L. Oppo, Phys. Rev. A 60, 2393 (1999). 
[8] R. Zambrini, M. Hoyuelos, A. Gatti, P. Colet, L. Lugiato, 

and M. San Miguel. Phys. Rev. A 62, 63801 (2000). 
[9] M. Santagiustina, P. Colet, M. San Miguel, and D. Wal- 
graef, Phys. Rev. Lett. 79, 3633 (1997). 

[10] M. Santagiustina, P. Colet, M. San Miguel, and D. Wal- 
graef, Opt. Lett. 23 1167 (1998). 

[11] R. J. Deissler, J. Stat. Phys. 40, 371 (1985); J. Stat. Phys. 
54, 1459 (1989); see also first work about convective sys- 
tems in plasma physics: R. J. Briggs, Electron- Stream 
Interaction with Plasmas, (M. I. T. Press, Cambridge, 
Massachussets, 1964). 

[12] M. Santagiustina, P. Colet, M. San Miguel, and D. Wal- 
graef, J. Opt. B 1, 191 (1999). 

[13] M. Santagiustina, P. Colet, M. San Miguel, and D. Wal- 
graef, Phys. Rev. E 58, 3843 (1998). 

[14] G. Iziis, M. Santagiustina, M. San Miguel, and P. Colet, 
J. Opt. Soc. Am. 16, 1592 (1999). 

[15] B. E. A. Saleh and M. C. Teich, Fundametals of Photon- 
ics, (Wiley, New York, 1991). 

[16] G-L. Oppo, M. Brambilla, L. A. Lugiato, Phys. Rev. A 
49, 2028 (1994). 

[17] H. Ward, M. N. Ouarazi, M. Taki and P. Glorieux, Eur. 
Phys. J. D 3, 275 (1998). 

[18] M. Taki, M. San Miguel and M. Santagiustina, Phys. 
Rev. E 61, 2133 (2000). 

[19] This spectrum is sharply narrowed on passing into the 
absolutely unstable regime |l|. 

[20] S. M. Barnett, P. Radmore, Methods in Theoretical 
Quantum Optics, (Oxford Science Publications, 1997). 

[21] C. W. Gardiner, P. Zoller, -2 enl. ed.- Quantum Noise, 
(Springer, Berlin, 2000). 



12 



[22] H. Carmichael, An Open Systems Approach to Quan- 
tum Optics, Lecture Notes in Physics, Vol. ml8, 
(Springer,Berlin, 1993). 

[23] R. F. Pawula, Phys. Rev. 162, 186 (1967). 

[24] A. Gilchrist, C. W. Gardiner, and P. D. Drummond, 
Phys. Rev. A 55, 3014 (1997). 

[25] W. H. Louisell, A. Yariv and A. E. Siegman, Phys. Rev. 
124, 1646 (1961). 

[26] S. Kumar and C. L. Mehta, Phys. Rev. A 21, 1573 (1980); 
ibid. 24, 1460 (1981) and references therein. See also S. 
M. Barnett and P. L. Knight, J. Opt. Soc. Am. B 2, 467 
(1985). 

[27] Numerical simulation are performed using a pseudo- 
spectral method as in We use the Gaussian random 
number generator proposed in: R. Toral, A. Chakrabarti, 
Comp. Phys. Comm., 74, 327 (1993). 

[28] The expression of supergaussian pump in all simulations 
is Eo(x) = e^a^/ 317 ) ; with x varying in the interval 
[—453; 453]. This produces a profile that is everywhere 
smooth but has a flat central plateau. 



[29] In presence of a flat pump, equations for ai(k,t) are in- 
variant under the change ai(k,t) — » e^ai^fc, i), so that 
it is equivalent to construct the W distribution consid- 
ering trajectories ai(k,t) or a'^k^t). For non flat pump 
this symmetry property is generally lost. In the pictures 
we represent the W distribution for the amplitudes ot\. 

[30] L. A. Wu, H. J. Kimble, J. L. Hall, and H. Wu, Phys. 
Rev. Lett. 57, 2520 (1986); S. Reynaud, Europhy. Lett. 
4, 427,(1987); A. Heidmann, et ai, Phys. Rev. Lett. 59, 
2555 (1987). 

[31] It is also possible to evaluate the squeezing outside the 
cavity. The output field arises from the interference of 
the intracavity field transmitted and the input field re- 
flected by the cavity input/output mirror (M.J. Collett 
and C.W. Gardiner, Phys. Rev. A 30, 1386 (1984)). The 
evaluation of the dynamics of the output field requires 
the simulation of coupled stochastic processes, as shown 




[32] J. Mertz, et ai, Opt. Lett. 16, 1234 (1991) 



13 



